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Abstract 

This paper focuses on cffectivity aspects of the Liiroth's theorem in differential 
fields. Let be a differential field of characteristic and F{u) be the field of differential 
rational functions generated by a single indeterminate u. Let be given non constant 
rational functions ui , . . . , w„ S J^{u) generating a subfield Q C F{u) . The differential 
Liiroth's theorem proved by Ritt in 1932 states that there exists v € Q such that 
Q = J-{v). Here we prove that the total order and degree of a generator v are bounded 
by miuj ord('yj) and {nd{e + 1) + l)^*^"*"^, respectively, where e := maxjOrd(wj) and 
d := maxj deg(vj). We also present a new probabilistic algorithm which computes the 
generator v with controlled complexity. 



1 Introduction 



In 1876, J. Liiroth in [IS] presented his famous result, currently known as the Liiroth's 
Theorem: if /c C L C k{u) is an extension of fields, where k{u) is the field of rational 
functions in one variable u, then L = k{v) for a suitable v (z L (see |27l §10.2] for a modern 
proof). In 1893 G. Castelunovo solved the same problem for rational function fields in two 
variables over an algebraically closed ground field. For three variables, Liiroth's problem 
has been solved negatively. 

In 1932 J.F. Ritt [22] addressed the differential version of this result: Let J- he a 
differential field of characteristic 0, u an indeterminate over T and T{u) the smallest 
field containing T , u and all its derivatives. Then, if G is a differential field such that 
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J- C G C J~{u), there is an element v ^ Q such that Q = J-{v). Such an element will be 
called a Liiroth generator of the extension J- <Z Q. 

In fact, Ritt considered the case of differential fields J- of meromorphic functions in an 
open set of the complex plane being Q a finite extension of J-. Later E. Kolchin in [TH] 
and |16j gave a new proof of this theorem for any differential field of characteristic and 
without the hypothesis of finiteness on Q. Contrary to the classical setting, the differential 
Liiroth problem fails does not hold the case of two variables (see [21] ) . 

The present paper deals with quantitative aspects of the Differential Liiroth's Theorem 
and the computation of a Liiroth generator u of a finite differential field extension. More 
precisely, let J-" be a differential field of characteristic 0, u differentially transcendental over 
J" and Q := T{Pi{u)/Qi{u), . . . , Pn{u)/Qn{u)), where Pj, Qj € T{u} are relatively prime 
differential polynomials of order at most e > 1 (i.e. at least one derivative of u occurs in Pj 
or Qj for some j) and total degree bounded by d such that Pj/Qj ^ for every 1 < j < n. 
First, we prove that any Liiroth generator v C Q can be written as the quotient of 
two relatively prime differential polynomials Mi (u) , M2 (u) € ^{u} with order bounded 
by min{ord(Pj/(5j); 1 < J < ''^j (see Proposition 1111 below) and total degree bounded by 
{nd{e + 1) + 1)2^+1 (see Section [O} • 

Secondly, we exhibit a probabilistic procedure, called Algorithm LurothGenerator 
(see Theorem [T6]) , which computes differential polynomials Mi (u) , M2 {u) such that its 
quotient is a Liiroth generator. If the input polynomials Pj, Qj, 1 < j < n, are given by a 
straight-line program of length L (for a definition and properties of straight-line programs 
see [2] and |13j). the number of arithmetic operations in performed by the algorithm is 
linear in L and polynomial in n,d,e and the degree of an algebraic variety V associated 
to the input polynomials (see Notation [T3] for the precise definition of V) . 

An algorithmic version of Ritt's proof of the differential Liiroth's Theorem is given 
in [8] but no quantitative questions on the order or the degree of the Liiroth generator 
are addressed and no complexity analysis of the algorithm is performed (for effectiveness 
considerations of the classical not differential version, we refer the interested reader to 

Our approach combines elements of Ritt's and Kolchin's proofs (mainly the intro- 
duction of the differential polynomial ideal related to the graph of the rational map 
u I— {Pj/Qj)i<j<n) with estimations concerning the order and the differentiation index 
of differential ideals developed in [23], [5] and [3- These estimations allow us to reduce 
the problem to a polynomial ring in finitely many variables. In this context, we are able 
to apply techniques of algorithmic resolution of algebraic polynomial systems a la Kro- 
necker as in [5] and [6] (see [12] and [26]) in order to obtain a probabilistic procedure with 
controlled complexity which computes the Liiroth generator. 

This paper is organized as follows. In Section 2 we introduce the notations, definitions 
and previous results from differential algebra (mainly concerning the order and the differ- 
entiation index) needed in the rest of the paper. In Section 3 we discuss some ingredients 
that appear in the classical proofs of Liiroth's Theorem by Ritt and Kolchin and that 
we will use in our arguments. In Section 4 and 5 we prove bounds for the order and 
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the degree, respectively, of any Liiroth generator. Section 6 is devoted to the algorithm 
and its complexity. Finally, in Section 7 we show how the algorithm works in two simple 
examples. 

2 Preliminaries 

In this section we introduce the notation used throughout the paper and recall some 
definitions and results from differential algebra. 

2.1 Basic definitions and notation 

A differential field {F, 6) is a field T equipped with a derivation 6 : J- ^ J-; for instance, 
J" = Q, M or C with 5 = 0, or J" = Q{t) with the usual derivation 6{t) = 1. 
Let {J-, 5) be a differential field of characteristic 0. 

The ring of differential polynomials in a indeterminates z := zi, . . . ,Za, which is de- 
noted by J-^{zi, . . . ,Za} or simply J-{z}, is defined as the commutative polynomial ring 
< j < a, p € No] (in infinitely many indeterminates), extending the derivation 
of J-" by letting 6{Zj^^) = Zj^~^^\ that is, zj*^ stands for the ith derivative of Zj (as custom- 
arily, the first derivatives are also denoted by zj). We write z^^^ := {z[^\ . . . ,z^'^} and 
zl^l := {z^^\ < i < p} for every p G Nq. 

The fraction field of J-{z} is a differential field, denoted by J-{z), with the derivation 
obtained by extending the derivation 5 to the quotients in the usual way. For g G ^{z}, 
the order of g with respect to Zj is ord(5, Zj) := max{i G Nq : -zj*^ appears in (7}, and the 
order of g is orA{g) := max{ord((7, Zj) : 1 < j < a}; this notion of order extends naturally 
to J-{z) by taking the maximum of the orders of the numerator and the denominator in a 
reduced representation of the rational fraction. 

Given a finite set of differential polynomials H = G J'{z}, we write [H] 

to denote the smallest differential ideal of J-{z} containing H (i.e. the smallest ideal 
containing the polynomials H and all their derivatives of arbitrary order). The minimum 
radical differential ideal of J-{z} containing H is denoted by {H}. For every i G N, we 
write i^W ■.= hf,...,hf and := H, H'^^\ . . . , H^^l 

A differential field extension QjT consists of two differential fields {J-,6jr) and {Q,6g) 
such that T Q and 5jr is the restriction to T of 5g. Given a subset T, C Q, denotes 
the minimal differential subfield of Q containing and S. 

An element G ^ is said to be differentially transcendental over if the family 
of its derivatives {^^^^ : p G No} is algebraically independent over otherwise, it is 
said to be differentially algebraic over T. A differential transcendence basis of Q/T is a 
minimal subset T, C Q such that the differential field extension Q/T{T) is differentially 
algebraic. All the differential transcendence bases of a differential field extension have 
the same cardinality (see |14t Ch. II, Sec. 9, Theorem 4]), which is called its differential 
transcendence degree. 
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2.2 Differential polynomials, ideals and manifolds 

Here we recall some definitions and properties concerning differential polynomials and 
their solutions. 

Let g € J^{z} = T{zi, . . . , Za}- The class of g for the order zi < Z2 < • • • < Za of the 
variables is defined to be the greatest j such that zj*^ appears in g for some i > if g ^ J-, 
and if 5 € If g is of class j > and of order p in zj, the separant of g, which will 

be denoted by Sg, is dg/dz^^ and the initial of g, denoted by Ig, is the coefficient of the 

(p) 

highest power of Zj in g. 

Given gi and g2 in J^{z}, g2 is said to be of higher rank in Zj than gi if either 

ovd{g2,Zj) > ord((jri,Zj) or ord(g2,-Zj) = ord(gi,Zj) = p and the degree of g2 in zj^^ is 

greater than the degree of 51 in Zj . Finally, g2 is said to be of higher rank than gi if g2 
is of higher class than gi or they are of the same class j > and g2 is of higher rank in Zj 
than gi. 

We will use some elementary facts of the well-known theory of characteristic sets. For 
the definitions and basic properties of rankings and characteristic sets, we refer the reader 
to O Ch. I, §8-10]. 

Let -ff be a (not necessarily finite) system of differential polynomials in J-{z}. The 
manifold of H is the set of all the zeros r/ € of for all possible differential extensions 

Every radical differential ideal {H} of J-{z} has a unique representation as a finite 
irredundant intersection of prime differential ideals, which are called the essential prime 
divisors of {H} (see [231 Ch. II, §16-17]). 

For a differential polynomial g in J-{z} of positive class and algebraically irreducible, 
there is only one essential prime divisor of {g} which does not contain Sg] the manifold of 
this prime differential ideal is called the general solution of g (see [231 Ch. II, §12-16]). 

2.3 Hilbert-Kolchin function and differentiation index 

Let ^ be a prime differential ideal of J-{z}. The differential dimension of denoted by 
diffdim(*p), is the differential transcendence degree of the extension ^ Frac(J^{z}/*P) 
(where Frac denotes the fraction field). The differential Hilbert-Kolchin function of ^ 
with respect to is the function Hr^ jr : No — >■ No defined as: 

Hij^ jr{i) ■= the (algebraic) transcendence degree of Frac(-F[2;W]/(*p n -F[zf*^])) over T. 

For i ^ 0, this function equals the linear function diffdim(*!p)(z + 1) + ord(*P), where 
ord(*p) € No is an invariant called the order of '^J ( |14l Ch. II, Sec. 12, Theorem 6]). The 
minimum i from which this equality holds is the Hilbert-Kolchin regularity of 

Let -F be a finite set of differential polynomials contained in *p of order bounded by a 
non- negative integer e. Throughout the paper we assume that e > 1, in other words, all 
the systems we consider are actually differential but no purely algebraic. 

Definition 1 The set F is quasi-regular at *p if, for every € No, the Jacobian matrix 
of the polynomials F,F,..., F^'^^ with respect to the variables has full row rank over 

the fraction field of J^{z}/^. 
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A fundamental invariant associated to ordinary differential algebraic equation systems 
is the differentiation index. There are several definitions of this notion (see [1] and the 
references given there) , but in every case it represents a measure of the implicitness of the 
given system. Here we will use the following definition, introduced in [U Section 3] , in the 
context of quasi-regular systems (i.e. finite sets of differential polynomials) with respect 
to a fixed prime differential ideal 

Definition 2 The *p— differentiation index a of a quasi-regular system F of polynomials 
in J-{z} of order at most e is 

a := min{A: G Nq : (F, F, . . . , F^''^)^^ n T\z^%^ = [F]<^ n T[z^%^, 

where := *P H (i.e. the contraction of the prime ideal and -F[z['^l]fp^ denotes 

the localized ring at the prime ideal • 

3 Differential Liiroth's Theorem 

In |22l Chapter VIII] (see also [14j and [23]), the classical Liiroth's Theorem for transcen- 
dental field extensions is generalized to the differential algebra framework: 

Theorem 3 (Differential Liiroth's Theorem) Let be a differential field of charac- 
teristic and let u be differentially transcendental over T . Let Q be a differential field 
such that F <Z Q C J^{u) . Then, there is an element v £ Q such that Q = F{v) . 

Our goal is the following: for n > 1, let be given differential polynomials Pi, . . . , P„, 
Qi> • • • > Qn £ •F{u}, with Pj/Qj ^ F and Pj,Qj relatively prime polynomials for every 
1 < J < "R, and denote 

g := F{Pi{u)/Qi{u), . . . , Pn{u)/Qn{u)), 

which is a subfield of T{u). We want to compute a Liiroth generator of Q/J-, that is, a 
pair of differential polynomials P,Q £ -F{n} such that Q ^ and Q = F{P{u)/Q{u)). 
We are also interested in the study of a priori upper bounds for the orders and degrees of 
both polynomials P and Q. 

We start by discussing some ingredients which appear in the classical proofs of Theorem 
[3] (see [231 US]) and that we also consider in our approach. 

Following [231 II. §39], let y be a new differential indeterminate and consider the 
differential prime ideal of all differential polynomials in Q{y} vanishing at u: 

S = {A G g{y} such that A{u) = 0}. 

Lemma 4 The manifold of S is the general solution of an irreducible differential polyno- 
mial B G Q{y}. More precisely, B is a differential polynomial in S with the lowest rank 
in y. 
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Proof. Let i? € S be a differential polynomial with the lowest rank in y. Note that 
B G G{y} is algebraically irreducible, since E is prime. We denote the order of i? by fc 
and the separant of B by Sb = dB/dy^^\ Consider the differential ideal 

^i{B):={AeG{y}\SBA = Q mod {5}}. 

As shown in j231 II. §12], the ideal Tii{B) is prime; moreover, we have that A G Tii{B) if 
and only if S'^A = mod [B] for some a € Nq and, in particular, if A G Hi[B) is of order 
at most k = ord(i?), then ^ is a multiple of B (see [23|. II. §13]). Furthermore, Tii{B) is 
an essential prime divisor of {B} and, in the representation of {B} as an intersection of 
its essential prime divisors, it is the only prime which does not contain Sb ( |23l II. §15]). 

In order to prove the lemma, it suffices to show that S = Si(i?). 

Let A € Tii{B). Then, SbA S {B}. Taking into account that S is prime, {B} C S, 
and Sb ^ S, it follows that A G S. 

To see the other inclusion, consider a differential polynomial yl G S. By the minimality 
of B, we have that A is of rank at least the rank of B. Reducing A modulo we obtain 
a relation of the type 

Sp%A = R mod [S], 

where h,c € No, Ib is the initial of B and R \s a differential polynomial whose rank is 
lower than the rank of B. Since A and B lie in the differential ideal S, it follows that 
i? G S, and so, the minimality of B implies that R = 0. In particular, S'^Ij^A G [B] and, 
therefore, G T,i{B). Now, Ib ^ T,i{B) since, otherwise, it would be a differential 
polynomial in S with a rank lower than the rank of B (recall that Si(i?) C S); it follows 
that A G Si(S). ■ 

Multiplying the polynomial B G G{y} by a suitable denominator, we obtain a differ- 
ential polynomial C G J-{u,y} with no factor in J-{u}. The following result is proved in 
[231 II. §42]: 

Lemma 5 If P{u) and Q{u) ^ are two coefficients of C (regarded as a polynomial in 
J-{u}{y}) such that P{u)/Q[u) ^ J-, the polynomial 

D{u,y) := Q{u)P{y) - P{u)Q{y) 

is a multiple of C by a factor in T , and Q = T{P{u)/Q{u)). ■ 

Note that, by the definition of C, the ratio between two coefficients of C coincides with 
the ratio of the corresponding coefficients of B. 

Under our assumptions, consider the map of differential algebras defined by 

ll) : T{xi,...,Xn,u} T{Pi{u)/Qi{u),...,Pn{u)/Qniu),u} 

Xi Pi{u)/Qi{u) 

li !->■ U 

Let *P C J-{x, ti} be the kernel of the morphism ip] then, we have an isomorphism 

F{Pl{u)/Qi{u), . . . ,Pn{u)/Qn{u),u} ~ J'jxi,. . .,Xn,u}/^. 
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This imphes that *p is a prime differential ideal. In addition, this isomorphism gives an 
inclusion 

T{Pi{u)/Qi{u), . . .,Pn{u)/Qn{u)} ^ T{xi, . . . ,x„,n}/<p, 

and the inclusion induced from this map in the fraction fields leads to the original extension 
Q = T{Pi{u)/Qi{u), . . . ,Pn{u)/Qn{u)) ^ ^{u) since *P does not contain non trivial 
polynomials in J^{u}. 

Now, if A S G{y} is a differential polynomial in S, multiplying it by an adequate ele- 
ment in T{Pi{u)/Qi{u), . . . ,Pn{u)/Qn{u)}, we obtain a differential polynomial in 
T{Pi{u)/Qi{u), . . . , Pn{u)/Qniu)}{y}, with the same rank in y as A. Taking a repre- 
sentative in J^{xi, . . . , Xn} for each of its coefficients, we get a differential polynomial 
A € J^{xi, . . . ,Xn,u}, with the same rank in y as A, such that A € *p. Conversely, 
from a differential polynomial M G J-{xi, . . . ,Xn,u} such that M € we may obtain a 
differential polynomial 

M{y) := M{Pi{u)/Qi{u),...,Pn{u)/Qn{u),y) eg{y} (1) 

vanishing at u, with a rank in y no higher than that of M. Notice that M G J'(nP<=])[y[2e]], 
since ord(M) < e and ord(Pj /Qj) < e for every 1 < j < n. 

We conclude that if M € T{xi, . . . ,Xn,u} is a differential polynomial in ^ with the 
lowest rank in u, the associated differential polynomial M{y) is a multiple by a factor in 
^ of the minimal polynomial B of u over Q. Therefore, a Liiroth generator of Q/J- can be 
obtained as the ratio of any pair of coefficients of M € g{y} provided this ratio does not 
lie in T. Moreover: 

Proposition 6 Let M € J-{xi, . . . , n} be a differential polynomial in ^ with the lowest 
rank in u. Consider two generic points ui,U2 G Q^'^^-^. Let Mi and M2 he the polynomials 
obtained by specializing M as in ([ip and substituting n'^^l for ui and U2 respectively in the 
polynomial obtained. Then Mi{u)/M2{u) is a Liiroth generator oj QjT. 

Proof. By Lemma [Sj we have that 

M(y) = M(Pi(n)/Qi(n), . . . , P„(u)/Q„(n), y) = ^{Q{u)P{y) - P{u)Q{y)) 

for some j £ G, and P,Q are such that Q = T{P{u)/Q{u)). Then, by means of two 
suitable specializations ui,U2 of the variables n'^^l (so that Qi, . . . ,Qn,Q and 7 do not 
vanish and P{ui)/Q{ui) 7^ Piu2)/Q{u2)), we obtain polynomials of the form 

Mi(y) = aiP{y) - PiQ{y) and M2{y) = a2P{y) - hQ{y). 

where ai , 02 ,/?!, /32 G and ai/32— a2/3i / 0. The proposition follows since J-'(Mi(u)/Af2(u)) = 
J^{P{u)/Q{u))=G. m 

For every j, 1 < j < n, denote Fj := Qj{u)xj — Pj{u) G J^{x,u}, and let F := 
Fi, . . . ,Fn. Then we have: 

Lemma 7 The ideal ^ is the (unique) minimal differential prime ideal of [F] which does 
not contain the product Qi . . . Qn- Moreover, ^ = [F] : (Qi . . . Qn)°° ■ 
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Proof. From the definitions of *P and F it is clear that [F] C ^ and Qi{u) . . . Qn{u) ^ ^■ 
Moreover, since is a characteristic set for the order in J-'{x, u} given by n < xi < • • • < x„ 
and *P n T{u} = (0), we conclude that for any polynomial if G *p there exists € No 
such that {Qi{u) . . . Qn{u))^ H € [F] (observe that Qj is the initial and the separant of 
the polynomial Fj for every j). The proposition follows. ■ 

Moreover, we have the following property that we will use in the sequel (recall Defini- 
tion [I]): 

Lemma 8 The system F is quasi-regular at *p. 

Proof. Let e be the maximum of the orders of the differential polynomials Fj, 1 < j < n. 
For every i € N, let Jj be the Jacobian matrix of the polynomials F^*"^! with respect to 
the variables (x, u)'*"^"'''^' . We have that 



and so, for every i € N, the minor of Jj corresponding to partial derivatives with respect 
to the variables x'*"^' is a scalar multiple of (Qi(n) . . . Qn{u)y, which is not zero modulo 
% ■ 

4 Bounding the orders 

Here we estimate the order in the variables x = xi, . . . ,Xn and n of a differential polynomial 
M(x, n) € ^ of minimal rank in u and the order of a Liiroth generator of the differential 
extension G/F. 

Let e := max{ord(Pj/(5j) : 1 < i < n}. Without loss of generality, we may assume 
that oid{Pi{u)/Qi{u)) > ••• > ovd{Pn{u)/Qn{u)). Consider the elimination order in 

F{xi , . . . ,Xn,u} with Xi < ■ ■ ■ < Xn < U. 

Taking into account that the fraction field of F{xi, . . . , Xn,u}/^ is isomorphic to 
F{u), we deduce that the differential dimension of *P equals 1. Since Pi{u)/Qi{u) ^ 
F, it is transcendental over F. Now, as the variable u^^^^^ appears in the derivative 
{Pi{u)/Qi{u))' but it does not appear in Pi{u)/Qi{u), it follows that {Pi{u) / Qi{u))' is 
(algebraically) transcendental over F{Pi{u)/Qi{u)). Continuing in the same way with the 
successive derivatives, we conclude that Pi{u)/Qi{u) is differentially transcendental over 
F. This implies that the differential ideal *P contains no differential polynomial involving 
only the variable Xi. 

Thus, a characteristic set of ^ for the considered elimination order is of the form 



Furthermore, Rn{xi, . . . differential polynomial in *P with a minimal rank in u, 

that is, we can take M{x,u) = Rn{x,u). Following [JU Lemma 19], we may assume this 




if h ^ j or h = j, k < I 
if /i = j, k > I 



i?l(xi,X2),i?2(xi,X2,X3), . . . ,i?„_i(xi, . . . , X„) , i?„ (xi , . . . ,X„, 



u). 
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characteristic set to be irreducible. Then, by [2A\ Theorem 24], we have that ord(i?i) < 
ord(*P) for every 1 < i < n; in particular, 

ord(M) < ord(q3). (2) 

The order of the differential prime ideal *P can be computed exactly: 

Proposition 9 The order of the differential ideal ^ equals e = max{oTd{Pj{u)/Qj{u)) : 
1 < i < n}. 

Proof. Lemma [7] states that the ideal *P is an essential prime divisor of [F], where 
F = Fi, . . . , Fn with Fj{x, u) := Qj{u)xj — Pj{u), 1 < j < n, and, as shown in Lemma [51 
the system F is quasi-regular at *p. Therefore, taking into account that e is the maximum 
of the orders of the polynomials in F, by [H Theorem 12], the regularity of the Hilbert- 
Kolchin function of ^ is at most e — 1. This implies that the order of ^ can be obtained 
from the value of this function at e — 1; more precisely, since the differential dimension of 
^ equals 1, we have that 

ord(<p) =trdeg^(-F[2;['^~iU''~^V(*Pn7-[x['=-iU[^-^l])) - e. 

In order to compute the transcendence degree involved in the above formula, we observe 
first that J^[x^''-^\u^''-^^]/{^ D J'[x^''-^\u^^-'^^]) ~ J^[{Pj /Qj)^''-^\u^^-^^. 

It is clear that the variables ul*^"^! are algebraically independent in this ring. Then, 
if L = T{u^^~^^), the order of the ideal *p coincides with the transcendence degree of 
L{{Pj/Qj)^^~^\j = l,...,n) over L. Without loss of generality, we may assume that 
ord(Pi/(5i) = e. Since the variable n^"^-* appears in Pi/Qi, we have that Liv!-^^) / L{Pi / Qi) 
is algebraic. Similarly, since n^*^"*"^) appears in (Pi/Qi)', it follows that the extension 
L(?/(<=),n(<^+^))/L((Pi/Qi),(Pi/Qi)') is algebraic. Proceeding in the same way with the 
successive derivatives of Pi/Qi, we conclude that L{u^^\ is algebraic 

over L((Pi/Qi)['^-i]). 

Since oxd{Pj/Qj) < e for j = 1, . . . , n, we have that L{{Pj / Q j = 1, . . . ,n) C 
L{u^^\u^^~^^\ . . . , n^^*^"^)) and, by the arguments in the previous paragraph, this extension 
is algebraic. Therefore, 

trdeg^L((P,/Q,)[^^, j = 1, . . . , n) = trdegiL(n(<^), . . . , = e. 



Then, by inequality ^ we conclude: 

Corollary 10 There is a differential polynomial M € *P C J-'{xi, . . . , x^, with the 
lowest rank in u such that ord(M) < e. 

Combining this corollary with the construction in Proposition El it is not too difficult 
to obtain the upper bound e for the order of the Liiroth generator Mi{u) / M2{u) of Q/J-. 
However, the following stronger result holds: 
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Proposition 11 Under the previous assumptions and notation, any element v £ Q such 
that Q = J^{v) satisfies ord(u) < mm{ord(Pj/(5j) : 1 < i < n}. 

Proof. Let u € ^ be such that Q = J-{v) and e := ord(u). 

For j = l,...,n, let Vj = Pj{u)/Qj{u). By assumption, vj ^ let Cj := OTd{vj). 
Let T be a new differential indeterminate over J^. Since Vj £ Q = y~{v), there exists 
Qj £ F{T) such that Vj = Qj{v). Let Nj = ord(ej). 

Assuming e > Cj, we have that Nj + e > ej and, therefore, 

du(^j+e) Z^^rw^ ^auW+^) 5t(^^)^ 

i>0 

Since A'j- is the order of @j, it follows that ^n^) 7^ 0, and as v is differentially transcen- 
dental over T, we have that -^^w^iv) ^ 0. We conclude that '^^1 = 0, which leads to 
a contradiction. ■ 

Note that the proof of the above proposition shows that all possible Liiroth generators 
V have the same order. In fact, two arbitrary generators are related by an homographic 
map with coefficients in (see for instance [15[ §1], \1'6\ Ch. II, §44]). 

5 Bounding the total degree 

5.1 Reduction to algebraic polynomial ideals 

As stated in Proposition[6l the Liiroth generator of Q /J- is closely related with a polynomial 
M{xi, . . . ,Xn,u) G *P with the lowest rank in u. 

By Corollarv llOl such a polynomial M can be found in the algebraic ideal ^fl ^-'[xt'^] , nt^l] 
of the polynomial ring u^'^l]. The following result will enable us to work with a 

finitely generated ideal given by known generators; the key point is the estimation of the 
^^-differentiation index of the system F := Fi, . . . , Fn (see Definition [2]) : 

Lemma 12 The ''^-differentiation index of F equals e. In particular, if := '^p n 
J^[x[^l,u[^l], we have that [F]^ n F[x^''\u^''%^ = {F, F , . . . , F^^^)^^ n J^ix^^'^u^^^. 

Proof. For every G N, let 3fc be the Jacobian submatrix of the polynomials F, . . . , F^^~^^ 
with respect to the variables {x,u)^^\ . . . , {x,u)^^^^~^\ The ^^-differentiation index of F 
can be obtained as the minimum k such that rank(^fc+i) — rank(5fc) = n holds, where the 
ranks are computed over the fraction field of F{x,u}/^ (see [H Section 3.1]). 

Now, since the order of the polynomials F in the variables x is zero, no derivative x^^^ 
with / > e appears in F, F , . . . , F^^~^\ This implies that the columns of the Jacobian 
submatrices of these systems corresponding to partial derivatives with respect to x^''\ 
I = e, . . . ,e + k—l, are null. On the other hand, as e is the order of the system F, we may 
suppose that the variable u^"^^ appears in the polynomial Fi and so, dFi/du^^^ ^ 0. Thus, 
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dF^^ /du'^^^ 7^ for /i — i = e and dF^^^ /du^^^ = for /i — i > e; that is, the matrices 5^, 



k = 1, . . . , e, are block, lower triangular matrices of the form 



/ ••• * 

••• * ••• * 



0*0 



* 







* / 



where denotes a zero column vector and * a non-zero column vector. Then, rank(^Jfc) = k 
for k = 1, . . . , e. Moreover, 



and, by the diagonal structure of 



(Qi{ 







Q2{u) 



\ 



••• * 



\ 



Qn{u) I 



we see that rank(5e+i) = rank(5e) + n. 

It follows that the *P-differentiation index of the system equals e. 



Notation 13 We denote byWc A^*^"*"^)" x A^^+^ the afftne variety defined as the Zariski 
closure of the solution set of the polynomial system F = 0, F = 0, . . . , F^^^ = 0, ni<j<n Qj 
0, where F = Fi, . . . , Fn and Fj{x, w^^]) = Qj{u^'^^)xj — Pj(u^^'^) for every 1 < j < n. 

The algebraic ideal corresponding to the variety V is {F, F, . . . , F^'^^) : , where 
q := Y\i<j<nQj- This ideal is prime, since it is the kernel of the map 



[e] [2el 
l<j<n ' 



Then, V is an irreducible variety. Moreover, F,F, . . . , F^'^^ is a reduced complete intersec- 
tion in {q 7^ 0} and so, the dimension of V is 2e + 1. With the notations of Lemma [T2\ we 
have: 

Lemma 14 The following equality of ideals holds: {F, F, . . . , F^^^) : q°° D F[x^^\ u'^'^^ = 
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Proof. First note that (F, F, . . . , : C (F, F, . . . , F(^))fp^ since qi'^. Conversely, 
like (f,f,...,f(^))^^ n 1], which is a prime ideal (see [H Proposition 3]), 

there is a polynomial g € F[x[^l,?x[ ^1] such that g and gh G (F, F, . . . , F^^)); but 

5 i (F, F, . . . , F(^)) : since otherwise, g G (F, F, . . . , F(^)) C *P for some G N 
contradicting the fact that q and 5 ^ Now, Lemma [12] implies that 

(F, F, . . . , F(^))qj^ n F[x[^Ut^l] = [Fh n n''^']. 
Finally, since [F]i;pF{x, n}fp = *PfpJ^{x, ti}fp (see [H Proposition 3]), we conclude that 

[F]^ n F[x[^U''^] = % n Ffxl'^U''^] = *P n u^"'] = *Pe. 

Therefore, 

(F,F,...,F(^)):g°°n.F[xW,u''^]=^e. 



5.2 Degree bounds 

From Proposition [6l in order to estimate the degree of a Liiroth generator it suffices to 
know the degree of the minimal polynomial M. For this, we will interpret M as an 
eliminating polynomial for the algebraic variety V under a suitable linear projection. 

Consider a finite set of variables Bq ^ constructed in the following way. Let 

Z he a maximal subset of {x'^^l} which is algebraically independent over T as elements of 
the field F(V) of rational functions over V, and U a maximal subset of such that 

{Z, U} is algebraically independent over F. Take Bq := {Z, U}. 

Proposition 15 The set Bq C {xt*^', ut^'} induces a transcendence basis o/F(V) over T. 
Proof. First, note that the canonical morphism 

is injective (Lemmata 1121 and I14|) . In particular Bq is a transcendence basis of the fraction 
field of F[x[^] , n F[xW, uW]). 

On the other hand, from [H Theorem 12], the order upper bound e is also an upper 
bound for the Hilbert regularity of Hence 

trdeg^F[xW, n T[x^^\u^^^]) = (e + 1) + ord(<P) = 2e + 1, 

as ord(*p) = e (see Proposition [9|), and therefore the cardinality of Bq is 2e + 1. 
Since 2e + 1 = dim(V) = trdeg_;rF(V) the proposition follows. ■ 

Let kf) := min{A; : u^'^^ ^ Bq}. By Corollary [lOl we have that /cq ^ e. Consider the 
projection 

^ : V ^ A2^+2, 7r(xW,n'^'') = {Z,U,u^''°'>). (3) 
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Then, from Proposition 1151 the Zariski closure of 7r(V) is a hypersurface in A^*^"*"^ and a 
square-free polynomial defining this hypersurface is the minimal polynomial M(x, u) we 
are looking for. In particular we have the inequality 

degM < degV 

(see im Lemma 2]) and we can exhibit purely syntactic degree bounds of V in terms of 
the number n of given generators for G/J-', their maximum order e, and an upper bound 
d for the degrees of their numerators and denominators. 

First, since the variety V is an irreducible component of the algebraic set defined by the 
{e + l)n polynomials F,F,..., F^^^ of total degrees bounded by d+ 1 (here F = Fi, . . . , Fn 
and Fj{x,u) = Qj{u)xj — Pj{u)), Bezout's theorem (see for instance |1H Theorem 1]) 
implies that 

degV< (d + l)(^+i)". 

An analysis of the particular structure of the system leads to a different upper bound 
for deg(V), which is not exponential in n: taking into account that V is an irreducible 
variety of dimension 2e + 1, its degree is the number of points in its intersection with a 
generic linear variety of codimension 2e + 1, that is, 

degV = #(VnF(Li,...,L2e+i)), 
where, for every 1 < i < 2e + 1, Lj is a generic affine linear form in the variables x^'^\u^'^^\ 

L,(xW,u''^l)= ^^ikxf+ b,ku^'^ + c,. (4) 

i<j<n 0<k<2e 

0<k<e 

For every 1 < i < n, the equation Fj(xj,u) = implies that, generic points of V 
satisfy Xj = Pj{u)/Qj{u). Proceeding inductively, it follows easily that, generically 

for every 1 < j < n, < /c < e. Substituting these formulae into (jH and clearing 
denominators, we deduce that the degree of Y equals the number of common solutions of 
the system defined by the 2e + 1 polynomials 



i<i<« 

i<j<n h^j 0<k<2e l<j<n l<j<n 

for generic coefficients aijk,bik,Ci and the inequality q ^ 0. From the upper bounds for 
the degrees of the polynomials Qj,Rjk, it follows that Lj is a polynomial of total degree 
bounded by nd{e + 1) + 1 for every 1 < i < 2e + l. Therefore, Bezout's bound implies that 



degV< {nd{e + + 



2e+l 
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Recalling that a Liiroth generator v of Q/F can be obtained as the quotient of two 
specializations of the polynomial M (see Proposition [6]) and that two arbitrary generators 
are related by an homographic map with coefficients in T (see |15l §1], \12>\ Ch. II, §44]), 
we conclude that the previous upper bounds for degV are also upper bounds for the 
degrees of the numerator and the denominator of any Liiroth generator of Q /J-. 

6 The algorithm 

This section is devoted to presenting an algorithm for the computation of a Liiroth gen- 
erator. The algorithm is based on the computation of the eliminating polynomial M 
associated to the projection vr defined in ([3|) over the positive-dimensional irreducible 
variety V. 

Our algorithm relies on the computation of a parametric geometric resolution of the 
system defining V as its main elimination step. Before proceeding, we briefly introduce 
this notion in our setting (for the general definition see Section 2.1]). 

Let B C {x''^!, lit^*^]} be a transcendence basis of T(y) over then, F{B) ^ T(V) 
is a finite algebraic extension. A parametric geometric resolution of V with respect to B 
consists in: 

• a linear form i in {rc^, tit^^^l} \ B which is a primitive element of J-{B) ^ -^(V), 

• the monic minimal polynomial m£ € J-{B)[T] of i, 

• parametrizations of the algebraic variables in terms of i, namely, for every variable 
w € {xW,u[^'^l} \ -B, a polynomial Vw € T{B)[T] with degV^ < degrrig such that 
^{£)w = K,(^) in -F(V). 

The data structure we use to represent multivariate polynomials is the straight-line 
program encoding. Roughly speaking, a straight-line program encoding a polynomial is 
a program which enables one to evaluate it at any given point. Each of the instructions 
in this program is an addition, subtraction or multiplication between two pre-calculated 
polynomials, or an addition or multiplication by a scalar. The number of instructions in 
the program is called the length of the straight-line program. For the precise definitions 
and basic properties we refer the reader to [2] (see also [l3]). 



Algorithm LurothGenerator 

INPUT: A positive integer e, a straight-line program of length L encoding polynomials 
Pi,Qi, . . . ,Pn,Qn G J'[n[''l] such that Q = F{Pi{u)/Qi{u), . . . ,Pn{u)/Qn{u)), and an 
upper bound d for the total degrees of the polynomials Pj, Qj. 

OUTPUT: An element v = Mi{u)/M2{u) G Q such that G = F{v). 

PROCEDURE: 

1. For J = 1, . . . , n, set Fj{x, u) := Qj{u)xj — Pj{u). 
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2. Compute an sip encoding F,F,..., F^^'^ and q := HKjXn Qii where F := Fi, . . . , 

3. Compute a transcendence basis q over 

with Z C x'*^] maximal and U C ut^'. 

4. Compute a parametric geometric resolution m, (Ko)t„g{a;W u[2<=l}\{z c/} ^ ■^(.^o)['^] of 
V with respect to the parameters Z, U. 

5. Set ko := min{/c : u^^^ ^ C/} and compute the minimal polynomial M{Z, U,u^^^^) of 
it(^o) modulo . . . : 

6. Prom two random specializations of M, obtain polynomials Mi(ti),M2(n) G J-'[n['^]] 
such that = F{Mi{u)/M2{u)). 



We now explain how the computations in the different steps of the algorithm are 
performed. 

In Step [H we use the input straight-line program in order to obtain a straight-line 
program of length L+^n encoding Fi, . . . , F„, q by simply computing, for each j = 1, . . . , n, 
a multiplication and a subtraction using the and previous results, 

and the product of the n polynomials Qi, • • • , Qn- 

From this straight-line program, in Step [21 we can obtain a straight-line program of 
length 0{e^n + e^L) encoding all the successive derivatives F,F, . . . , F^^^ proceeding as in 
O Lemma 21] (see also [191 Section 5.2]). 

The computation of the transcendence basis Bq = {Z,U} involved in Step [3] relies on 
the well-known Jacobian criterion from commutative algebra [171 Chapter VI, Section 1, 
Theorem 1.15]; more precisely, we apply [Jl Lemma 19]. Since J-[x^^\u^'^^'^]/{F, F, . . . , F^'^^) : 
qoo ^ F[{Pi{u)/Qi{u))^''\. . . , (P„(ti)/Q„(n))['^],ti[2<^]], it follows that ah rank computations 
can be done in J^(n[^'^l). Therefore, we may proceed as in the proof of [51 Theorem 27]. 
The resulting subroutine is polynomial in e, n, and linear in L. 

Step m is achieved by means of a slightly modified version of the algorithm in |261 
Section 4] within a complexity polynomial in n, e, d, deg V and linear in L. The algorithm 
proceeds in three main steps: 

(a) Initial estimation: Given a point (C^^) ^ Q^^"*"^, compute a geometric resolution of 
the simple roots of the system F|((^ = 0, = 0, . . . , F^^^ |((^ = 0, 7^ 
obtained by specializing (Z, U) = ((", i^)- 

(b) Approximation: Starting from the solution of the specialized system, by means of 
a Newton-Hensel type procedure approximate the coefficients of the corresponding 
parametric resolution in a ring of formal power series. 

(c) Reconstruction: Recover the coefficients of the parametric geometric resolution from 
their power series expansion at {C,i^)- 

In fact, our modification of the procedure in [261 Section 4] consists in the way we ma- 
nipulate the objects in intermediate computations: we use straight-line programs to deal 
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with truncated power series. In particular, we modify step (c) of the algorithm, which is 
achieved by a multivariate Pade approximation procedure ( \2()\ Section 4.3.1]), in order to 
work with the straight-line program encoding of polynomials, by replacing the Euclidean 
extended algorithm with subresultant computations (see [3 Section 5.9, Corollary 6.49]) 
as explained in the proof of O Proposition 46]. 

Once the parametric geometric resolution of Y is obtained, the computation of the 
minimal polynomial in Step [5] is straightforward (see, for instance, |12j ) and can be done 
within the same complexity bounds. 

We point out that the computation of the polynomial M could be performed, alterna- 
tively, from the polynomials F,F,..., F^'^\q by a procedure based on an extension of the 
method presented in |12j (which assumes the projection to be a finite morphism) to the 
case of a dominant map by means of the techniques described in [2^ (this approach can 
be found, for instance, in [5l Proposition 46]). 

Finally, in Step [6] we specialize the variables x and its derivatives in M[x,u) into the 
values obtained by choosing at random two sets ui,U2 of values for ut^^l and solving the 
linear system F{x, Ui) = 0, F{x, x,Ui) = 0, . . . , F^^\x, x, . . . , x^'^\ui) = in the unknowns 
x^^\ In this way, we obtain polynomials Mi{y), M2{y) providing a Liiroth generator v := 
Mi{u)/M2{u) of G/J~- Note that the matrices of the linear systems to be solved are the 
specializations at ui,U2 of the Jacobian matrices of F,F, . . . jF^^') with respect to the 
variables x^'^'^ with an additional column containing F{0,Ui), . . . , F^^\o,Ui) and so, they 
can be computed from a straight-line program encoding the polynomials in polynomial 
time. 

Therefore, we have the following complexity result: 

Theorem 16 Let T be a differential field of characteristic and u differentially tran- 
scendental over F. Let Q = F{Pi{u)/Qi{u), . . . , Pn{u)/Qn{u)) , where Pj,Qj € F{u} are 
relatively prime differential polynomials of order at most e and total degree bounded by d 
such that Pj/Qj ^ F for every 1 < j < n. Algorithm LurothGenerator computes differen- 
tial polynomials Mi(u), M2{u) of order and total degree bounded by min{ord(Pj/(5j); 1 < 
j < n} and degV respectively, such that Q = F {Mi{u) / M2{u)) , where V is the affine 
variety introduced in NotationM^ 

From a straight-line program of length L encoding the polynomials Pj,Qj, 1 < j < 
n, the algorithm produces a straight-line program encoding Mi,M2. The length of the 
output straight-line program and the number of arithmetic operations in F performed by 
the algorithm are linear in L and polynomial in n, d, e and deg V. 

7 Examples 

As before, let be a differential field of characteristic and u a differentially transcen- 
dental element over F. 

Example 1. Let Q = F{u/u, u-\- ii). In this case, we have: 

• e = 1, n = 2, 

• Fi = iixi — u, F2 = X2 — u — it,, q = u, 
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As the ideal {Fi, ^2,^1,^2) = (^2:1 — u, X2 — u — u, ilxi +uxi — u,X2—u — u) is prime and 
does not contain ii, we have that (^1,^2,^1,^2) : 9°° = {Fi, F2, Fi, F2) and, therefore, 

V = V{Fi, F2, Fi, -F2) = V{uxi — u, X2 — u — ii, ilxi + iixi — u,X2 — u — u). 

The dimension of V equals 3. We obtain a suitable transcendence basis of the fraction 
field of F[xi, X2,xi,X2,u, u, u]/{Fi, F2, Fi, -F2) over by studying the Jacobian matrix of 
the polynomials: 







(ii 











— 1 Xl 




(ii 











— 1 u/ii 


\ 


5(Fi 


F2,Fi,F2) _ 





1 








-1 -1 








1 








-1 -1 





d{xi,X2 


Xi,X2,U,U,u) 


ii 





ii 





Xl - 1 


Xl 


~ u 





ii 





—uujii} 


u/ii 














1 


-1 


-1/ 


Vo 








1 


-1 


-1/ 



Since the matrix obtained by removing the first three columns of this last matrix has full 
rank in F{u,u,ii), we deduce that {xi,X2,xi} is a transcendence basis as required. 

Then, ko = and so, we look for a polynomial M(xi, X2, xi, n) G {Fi, F2, Fi, F2). We 
compute this polynomial by an elimination procedure: 



M{xi,X2,Xi,u) = {Xi + l)u - X\X2. 

Finally, specializing {u,ii,u) into ui = (1,1,0) and U2 = (0,1,0), we compute specializa- 
tion points (1, 2, 1) and (0, 1, 1) respectively for (xi, X2, xi); hence, we obtain the following 
Liiroth generator for Q/J-': 

_ Mi(m) _2u-2 
^ ~ M2{u) ~ u 

In the previous example, the minimal polynomial M lies in the polynomial ideal 
(-Fi,-F2)i that is, no differentiation of the equations is needed in order to compute it 
and, consequently, a Liiroth generator can be obtained as an algebraic rational function 
of the given generators. However, this is not always the the following example 

shows. 

Example 2. Let Q = F{ii,u + ii). We have: 

• e = 2, n = 2, 

• Fi = Xl — u, F2 = X2 — u — u^'^\ q = 1, 
Following our algorithmic procedure, we consider the ideal 

(Fi , F2 , Fi , F2 , Ff ^ , Ff ) ) = (xi -it, X2 -n-n(2) , ii -u^^) , ±2 -n-n^^) , xf'> -u^^\ xf -n^^) -u^^) ) . 

This is a prime ideal of J-'[x[^], u^]. The variety V is the zero-set of this ideal and it has 
dimension 5. A transcendence basis of F(V) over T including a maximal subset of x^ is 
{xi, X2, Xl, X2 5 2^2^''} and the minimal polynomial of u over F(xi, X2, xi, X2 5 2:2^'') is 

M = — X2 + Xl. 

Two specializations of this polynomial lead us to a Liiroth generator of Q jJ- of the form 
v = ^. We conclude that Q = J'iu). 
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